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Abstract 

We discuss the nuclear interactions of neutrinos versus those of antineutrinos, a relevant compar¬ 
ison for CP violation experiments in the neutrino sector. We consider the MiniBooNE quasielastic¬ 
like double differential neutrinos and antineutrinos cross sections which are flux dependent and 
hence specific to the MiniBooNE set-up. We combine them introducing their sum (v + v) and their 
difference (is — is). We show that the last combination can bring a general information, which can 
be exploited in other experiments, on the nuclear matrix elements of the axial vector interference 
term. Our theoretical model reproduces well the two cross sections combinations. This confirms 
the need for a sizeable multinucleon component in particular in the interference term. 

PACS numbers: 25.30.Pt, 13.15.+g, 24.10.Cn 
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I. INTRODUCTION 


One of the challenging goals of neutrino experiments is the detection of CP violation in 
the neutrino sector. A convincing evidence would be the detection of an asymmetry between 
the oscillation rates of muon neutrinos and antineutrinos into electron neutrinos since, in the 
absence of CP violation, these rates are the same. For this test one needs a muonic neutrino 
beam and an antineutrino one. The electrons (positrons) produced in a far detector by the 
charged current interaction of the electron neutrinos (antineutrinos) are the signature for 
the —y u e (i? M —> v e ) oscillation process, once the direct v e (v e ) background is eliminated. 
Several obstacles can stand on the way of the detection of CP violation through the u, v 
asymmetry. One is that the interactions of neutrinos and antineutrinos with any nucleus 
are not identical but they differ by the sign of the axial vector interference term, creating an 
asymmetry unrelated to CP violation and which has to be fully mastered. This is not a trivial 
task due to the complexity of the nuclear dynamics. It reflects in the neutrino interactions 
and may obscure the message that one wants to extract on the oscillation mechanism. 
One example concerns the role of the multinucleon emission process which in a Cherenkov 
detector is misidentihed for a quasielastic one [1, 2], This misinterpretation produces an 
apparent increase of the neutrino quasiclastic cross section, at the origin of the so called 
axial mass anomaly detected in the MiniBooNE experiments [ 3 J] as now widely accepted 


El 4- 


15]. The detection of CP violation which involves a comparison between neutrinos 


and antineutrinos events needs an even more detailed understanding of the multinucleon 
processes since it concerns the difference between neutrino and antineutrino cross sections. 
This is not a priori obvious and it is the object of the present article. 


II. ANALYSIS AND RESULTS 

In order to illustrate the respective roles of the multinucleon components in the v and v 
cross sections, as we have introduced in a previous work [ 4 ], we start by giving below the 
following simplified expression (we remind however that for the actual evaluation we use a 
more complete one) for the double differential neutrino or antineutrino cross sections on a 
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nuclear target such as 12 C 


d 2 a 



d cos Qduj 


( 1 ) 


where the plus (minus) sign applies to neutrinos (antineutrinos). In this Eq. (ITj) Gp is the 
weak coupling constant, 6 C the Cabibbo angle, Ge, Gm and Ga are the nucleon electric, 
magnetic and axial form factors, E v and Ei the initial and final lepton energies, ki the 
modulus of the final lepton momentum, c o and q the energy and the momentum transferred 
to the nucleus and 6 is the lepton scattering angle. The cross section on the nuclear target 
is expressed here in terms of the nuclear responses R to probes with various couplings to the 
nucleon, isovector (index r) or isovector with isospin and spin coupling (index err). For the 
last responses the isovector spin coupling can be spin transverse, crr(T), or spin longitudinal, 
ctt(L). The responses are q and cu dependent. The last term of Eq. ([T]) which changes sign 
between v and v is the axial vector interference term, the basic asymmetry which follows 
from the weak interaction theory. It is expressed here in terms of isospin spin-transverse 
nuclear response. Its evaluation for complex nuclei is not trivial due to the presence of 
many-body effects. This is why it is important to obtain an experimental information on 
this term, one object of the present work. 

Now, for our description of the multinucleon component of the responses, we followed 
the experimental indications provided by the electron scattering data. In the transverse, 
or magnetic, data the dip between the quasielastic and the Delta part of the response 
is filled, which we interpreted as an indication for the presence of two nucleon emission 
(l6 |. In the charge response instead this component is absent. With these indications, for 
neutrinos we have introduced the multinucleon component only in the spin isospin response, 
the pure isovector one keeping its quasielastic character. This difference is at the origin 
of the statement in our work of Ref. Q] that the multinucleon effect should be relatively 
less pronounced for antineutrinos than for neutrinos, the reason being that the isovector 
response which is free from multinucleon effect plays a larger relative role for antineutrinos 
due to the negative sign of the interference term. This argument obviously holds only in the 
kinematical regions where the influence of the isovector component is significant. 
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We have previously tested our model on the MiniBooNE data for the differential cross 

8] and antineutrinos 13]. The good fit of the 


sections 


17], independently for neutrinos 


data provided by our model did not contradict our statement on the respective roles of 
np-nli. However in these works, the test was performed separately for the neutrinos and the 
antineutrinos and we did not specifically address the detailed comparison between the two. 
As this comparison is essential for the CP violation detection we want to investigate in the 
present work this question in more details. For Cherenkov detectors, such as the MiniBooNE 
or the Super-Kamiokande ones, are we able to describe quantitatively the difference between 
neutrinos and antineutrinos nuclear interactions? The wealth of experimental data which has 
accumulated in the last years by the MiniBooNE experiment allows an experimental test for 
this comparison and we will explore the message which can be drawn from these data. The 
quantity best suited for this exploration is the measured double differential quasiclastic-like 
(i.e. which incorporates also multinucleon states) cross section ,, d ^ dF with respect to the 




muon energy E M and the muon emission an; 
problem of the neutrino energy 181], |l 1]. [19 

d 2 a f 


c 9, which is not affected by the reconstruction 


2 If . However it is flux dependent, with 


d cos 9dE„ 


d 2 c 
d cos 6du 


j=e v —Eu §{E v )dE u 


( 2 ) 


where oj is the energy transferred to the nucleus and <&(E U ) is the neutrino (or antineutrino) 


normalized flux. Neutrino experiments measure 


d?c 


d cos OdEn 


while nuclear physics evaluations 


calculate the quantity 


d?c 


d cos Odw ’ 
d?c 


which is the basic ingredient for all analysis on neutrino data. 
The flux dependence in dc ^ ddE is a priori an obstacle to extract a universal comparison 
between the two cross sections for neutrinos and antineutrinos, applicable to CP violation 
experiments. For each set of flux profiles the measured differential cross sections are different. 
In particular the v and D asymmetry for dc ^ ddE has two sources, one arises from the 
basic weak interaction theory, as given in Eq. (IT]) . The second one which arises from the 
flux asymmetry has no universal character and is specific for each experiment. Is there 
nevertheless something general and informative in the MiniBooNE data? It is the question 
that we will address in the following. 

Let us consider the following combinations: the sum, sum, and the difference, dif , of the 
double differential cross sections for neutrinos and antineutrinos with respect to the lepton 
emission angle, 9, and to the energy c o transferred to the nucleus 

d 2 a u d 2 a D 


sum(cos9,u) = 


dcos9du d cos 9du 


( 3 ) 


4 





















E v (GeV) 


FIG. 1: (color online) Normalized MiniBooNE v^ and fluxes derived from Refs. [3|] and 
respectively. Their half sum ( < h+) and half difference (<!>_) are also shown. 


Q 


while for the difference dif the sign plus is changed to a minus one. Notice from Eq. ([TJ) 
that the difference, dif, contains only one term, the axial vector interference one. It is this 
quantity which governs the difference in cross sections of v and v but it is not accessible 
directly from neutrino data since only derivatives with respect to the emitted lepton energy 
are measurable. We also introduce S and D , the corresponding flux integrated quantities, 
which are instead experimentally accessible in the MiniBooNE data, 


S(cos 9, E,f) = 


D (cos O.E^) = 


d 2 


cr v 


+ 


d 2 a p 


d cos 9dE u d cos OdE, 


d 2 a,. 


d 2 a D 


(4) 

(5) 

d cos 6dEf d cos 6dEf 

in which quantities such as dc f s g dE have been defined previously in Eq. (121) . The MiniBooNE 
experimental distribution of D(cos 9, E^) has been given by Grange and Katori 22] in a 
tridimensional plot. 

If the normalized neutrino and antineutrino flux profiles would be identical, ^ l/ {Ef) = 
^q(E v ), their common value could be factorized in the integrals over the neutrino energies 
implicitly contained in the above Eqs. (J4]) and (l5l) . In this case only the axial vector 
interference term would survive in the difference D, while the sum S would totally eliminate 
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this term. This would allow a direct experimental test of the interference part. However this 
is not quite the case as appears in Fig. [T] which compares the two normalized MiniBooNE 
fluxes. In view of this difference we express the two cross section combinations in terms 
of the average flux, $ + (E u ) = 1/2 [Q y (E u ) + $ P (E V )\ and of the flux difference, $_(_£?„) = 
1/2 [$ U (E U ) — <1 ) D (E U )\ as follows 

S(cosd,E^) = I dE u [sum(cos9,u})\ u=Ev -E^+(E v ) + dif (cos9, uj)\u =e „-e^-(E„)] (6) 

and 


D{[cos6, Eft) = J dE v [sum(cosd,U})\ u=Ev -E^-(K) + dif(cos9,U})\ u=Ev -E^+(E v )]. (7) 

We observe again that for identical fluxes, for which <&_(E V ) = 0, the quantity D probes only 
the quantity dif , i.e., the axial vector interference term, while the sum S is totally blind to 
it. Notice that the MiniBooNE flux difference <fr-(E u ) is negative for small E v values and 
becomes positive beyond an energy ~ .7 GeV. 

The question is now the following: in view of the fluxes difference, what remains in the 
MiniBooNE data of the purity of the difference D with respect to the axial vector interference 
term and accordingly of its elimination in the sum S'? In order to answer these questions 
we have evaluated the sum S and the difference D with the real MiniBooNE fluxes on 
the one hand and with only the mean flux $+(.£/,) on the other hand. These evaluations 
are performed with our theoretical model which has been described for instance in Refs. 
( 2 , 8|. It incorporates the multinucleon component of the cross section and the responses 
are treated in the random phase approximation, RPA. The resulting curves are displayed 
in Figs. \Z\ and 0 as functions of the muon angles, for various values of the muon kinetic 
energy. For what concerns the sum S, for all angles and in the full range of muon energies, 
the two sets of curves are very close, confirming the major role of the mean flux contribution 
and showing that there is very little sensitivity in the sum on the axial vector interference 
term. For the difference D the two sets of curves are quite close up to a muon kinetic 
energy of X/ ~ .6 to .7 GeV. Beyond this energy they progressively depart. The fact that 
the mean flux curve is below the genuine flux one reflects the positive sign of $_(£/,) for 
large E u values (large X/ values meaning also large E u values). Notice also that for the 
smallest 7/ values, for which the neutrino energy can be small, one can detect a change in 
the sign of the effect, which becomes slightly negative, reflecting the sign change of $_(£/,). 
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FIG. 2: (color online) Sum of the v and v flux folded double differential cross sections on carbon 
per active nucleon plotted as a function of the cos 0 for different values of the emitted muon kinetic 
energies. Continuous line: evaluation (RPA + np-nh) with the real v and v MiniBooNE fluxes; 
dot-dashed line: evaluation (RPA + np-nh) with the mean flux d>_|_. 


In spite of the (moderate) deviations at large muon kinetic energies, for all T p values the 
contribution from the mean flux remains dominant in the difference D. This is fortunate 
since it means that the axial vector interference term is experimentally accessible in the 
MiniBooNE data through the neutrino and antineutrino double differential cross sections 
difference. We can in particular test if our multinucleon component, which in our model is 
maximum in the axial vector interference term, i.e., in the difference D , is compatible, or 
not, with the data. While with the sum of the neutrino antineutrino cross sections instead 
we can explore its role in the remaining part of the cross section. In general when the v 
and v normalized flux difference is small the same conclusion applies: the difference in the 
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FIG. 3: (color online) Difference of the v and v flux folded double differential cross sections on 
carbon per active nucleon plotted as a function of the cos 9 for different values of the emitted muon 
kinetic energies. Continuous line: evaluation (RPA + np-nh) with the real u and v MiniBooNE 
fluxes; dot-dashed line: evaluation (RPA + np-nh) with the mean flux 4> + . 


measurable differential cross sections between neutrinos and antineutrinos is dominated by 


the axial-vector interference term. This is the case for T2K 


23 ] and also for the NuMI 24] 


beams, the ones used in the MINOS, MINERVA and NOz/A experiments. 

Having seen the message that the comparison with the mean flux curves carries i. e. that 
the axial vector interference term dominates the difference and has very little influence on 
the sum, in the following we calculate in our model the sum and the difference of the cross 
sections with the true neutrinos and antineutrinos fluxes, and not with the mean one, in 
order to avoid unnecessary errors. Our present results are then simply the sum and the 
difference of our previous theoretical results published in Ref. js] for neutrinos and in Ref. 
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FIG. 4: (color online) Sum S of the v and v MiniBooNE flux folded double differential cross sections 
on carbon per active nucleon plotted as a function of the cos 6 for different values of the emitted 
muon kinetic energies. Continuous line: our complete RPA evaluation including the multinucleon 
emission channel; dashed line: genuine quasielastic contribution calculated in RPA; dot-dashed 
line: quasielastic contribution in the bare case. The points are the combination of the MiniBooNE 
experimental results 3.1171. 


Q for antineutrinos. In Figs. [4] and [5] we display our calculated values of ,S'(cos 0, A /( ) and 
_D(cos$, E^) as a function of the muon emission angle, for various values of the muon kinetic 
energy, together with the experimental data points. Our predictions which incorporate the 
multinucleon component account quite well for the data for all angles and in the full range 
of muon energies, both for the sum and for the difference. Only for the smallest 7j t bin some 
small deviations can be observed. 

This test of the validity of our model for these combinations is important as it addresses 
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FIG. 5: (color online) The same of Fig. [I] but for the difference D. 


directly to the question of our understanding of the neutrinos versus antineutrinos interac¬ 
tions, in particular for what concerns the crucial and more debated role of the multinucleon 
component. Our predicted cross sections without the multinucleon part are also shown in 
Figs. [4] and 0 they definitely fail to account for the data both for the sum and the difference. 
The necessity of the introduction of a large multinucleon component, in particular in the 
difference i.e. in the axial vector interference term, is confirmed here by the curves of Fig. 
[5] Another check remains to be performed to fully confirm the necessity of the multinucleon 
piece, in particular in the axial vector interference term: in our description the collectivity 
of the quasielastic responses has been included in the form of the random phase approxi¬ 
mation treatment, RPA. It produces a suppression effect of the quasielastic part due to the 


repulsive character of the particle-hole force 


251 ]. For the spin-isospin transverse response 


it is the Ericson-Ericson-Lorentz-Lorenz quenching arising from the mixture of Delta-hole 
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FIG. 6: (color online) Sum and difference of the v and v flux folded differential cross sections 
d ^ s9 and jjjr- on carbon per active nucleon. Continuous line: complete RPA evaluation including 
the multinucleon emission channel with the real u and u MiniBooNE fluxes; dashed line: genuine 
quasielastic contribution calculated in RPA with the real v and u MiniBooNE fluxes; dot-dashed 
line: complete RPA evaluation including the multinucleon emission channel with the mean flux 
<f> + ; dotted line: contribution of the isovector response only. 


states into nucleon-hole ones 


26]. The question is then: if a good fit is obtained with the 


combined and opposite effects of RPA and of the multinucleon component, could it be that a 
similar good fit would be achieved by omitting both effects. Would the simplest quasielastic 
description also account for the data? In Figs. Q] and [21 the effect of RPA is suppressed 
in the quasielastic cross section, which indeed has some enhancement effect. In particular 
for the cross section difference this enhancement is moderate and not enough to account by 
itself, for the data. We can safely conclude that a large multinucleon component is needed 
to describe the data for the axial vector interference term which governs the cross sections 
difference. The same conclusion applies to the cross sections sum, i.e. to the remaining 
part of the interaction, it is also appreciably influenced by the multinucleon component. 
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FIG. 7: (color online) Sum and difference of the v and v flux folded Q 2 distributions on carbon per 
active nucleon. Continuous line: complete RPA evaluation including the multinucleon emission 
channel with the real v and v MiniBooNE fluxes; dashed line: genuine quasielastic contribution 
calculated in RPA with the real v and u MiniBooNE fluxes; thin dot-dashed line: quasielastic 
contribution calculated in the bare case with the real u and u MiniBooNE fluxes; thick dot-dashed 
line: complete RPA evaluation including the multinucleon emission channel with the mean flux 
$+• 


Our model for the neutrino nucleus interaction is able to describe both components. As an 
additional illustration of its success we report in Fig. [6] the single differential cross sections, 
with respect to the muon kinetic energy or to the muon emission angle. As previously we 
deal with the sum and the difference for neutrinos and antineutrinos calculated with the true 
and the averages fluxes. We can observe again that the sum shows practically no sensitivity 
to the flux difference, while the cross sections difference displays some a mild sensitivity, in 
particular in the forward direction or at large T ^ values. Our predictions which include the 
multinucleon component reproduce well the data. In our model the relative importance of 
the multinucleon term in the cross sections combinations depends on the role of the isovector 
response which is shown in Fig. [6l the smaller this role the larger the multinucleon contribu¬ 
tion. It is the largest in the cross section difference. Similarly the isovector response weight 
is larger in antineutrino cross sections than in neutrino ones, hence the smaller multinucleon 
contribution for antineutrinos. It is however not a large difference and the multinucleon 
influence remains important also for antineutrinos. 
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Finally for completeness we combine in Fig. [7] our previous evaluations of the neutrino 
and antineutrino Q 2 distributions, published in Refs. [8] and [13] respectively, to evaluate 
their sum and difference. We also display the result obtained with the averaged flux $ + . 
As previously pointed out the difference is more sensitive to the mean flux approximation 
which however gives the bulk of this difference. As a consequence the difference of the 
experimental MiniBooNE points 
of the axial vector interference term. This conclusion would apply as well to the MINERVA 
neutrino [27] and antineutrino 28] Q 2 distributions, due to the closeness of the neutrino and 


da 

dQ 2 


Ap ) is directly related to the Q 2 distribution 


antineutrino normalized fluxes. It will be the object of a future investigation. 

Notice that the present study is done for the muonic neutrinos of the MiniBooNE exper¬ 
iment, for which data are available while CP violating experiments through the asymmetry 
of the oscillations rates of muonic neutrinos and antineutrinos into electron ones involve 
the detection of electrons or positrons produced in a detector by the charged current in¬ 
teractions of these electron neutrinos. We have already addressed the question of electron 
neutrino cross sections, versus the muon neutrino ones 21]. The effect of the small change 
in kinematics due to the smaller electron mass will not affect the present conclusions: the 
electron neutrino nuclear interactions produce by themselves an important v e v e asymmetry. 
The present study indicates that the multinucleon role is essential in this problem and to 
what precision this asymmetry can be mastered. 


III. SUMMARY AND CONCLUSIONS 

In summary we have investigated two different combinations of the neutrino and antineu¬ 
trino MiniBooNE flux folded double differential cross sections on 12 C, their sum and their 
difference. They probe different pieces of the neutrino or antineutrino interactions with 
the nucleus. These quantities depend on the neutrino or antineutrino flux. In the case of 
identical neutrino and antineutrino normalized flux profiles, the difference provides a direct 
access to the axial vector interference term, while the sum totally eliminates it. As this is 
not actually the case for the MiniBooNE fluxes we have tested how much the flux difference 
influences the two combinations. We have shown that this influence is small both on the 
sum and on the difference of the cross sections. These combinations remain rather pure 
with respect to the axial vector interference term which is either dominant (difference) or 
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nearly absent (sum). This allows more specific tests of our theoretical model on the axial 
vector interference term, important for the CP violation data. Our model gives a good fit 
for the MiniBooNE data for the sum and the difference of the cross sections reproducing 
well the data in the full range of muon energy and emission angle. The introduction of 
the multinucleon component is necessary for a good fit, an important test for its presence 
in the axial vector interference term. The success of our description indicates that we can 
reach a good understanding of the nuclear effects in neutrino interactions, also for what 
concerns the comparison between neutrinos and antineutrinos cross sections. We have con¬ 
centrated in this work on the interactions of muonic neutrinos where a complete set of data 
is available. Our predictions can easily be extended to electron neutrinos, relevant for CP 
violation data. The nuclear cross section difference for neutrinos and antineutrinos stands 
as a potential obstacle in the interpretation of experiments aimed at the measurement of 
the CP violation angle, 6. The present analysis, performed on the MiniBooNE data, shows 
the importance, for Cherenkov detectors, of the inclusion of the multinucleon contribution 
to the quasielastic-like cross section for mastering this difference. 
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